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We study the reach of the Large Hadron Collider with 1 fb -1 of data at yfs = 7 
TeV for several classes of supersymmetric models with compressed mass spectra, 
using jets and missing transverse energy cuts like those employed by ATLAS for 
Summer 2011 data. In the limit of extreme compression, the best limits come from 
signal regions that do not require more than 2 or 3 jets and that remove backgrounds 
by requiring more missing energy rather than higher effective mass. 
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I. INTRODUCTION 

The Large Hadron Collider (LHC) is now testing the proposal that super symmetry [l| (SUSY) is 
the solution to the hierarchy problem associated with the electroweak scale. At this writing, there 
have been no hints of supersymmetry, defying expectations based on the sensitivity of the Higgs 
potential to superpartner masses in many models, including the popular "mSUGRA" (minimal 
supergravity) scenario. It is possible that the failure of SUSY to appear is simply due to the up 
and down squarks being very heavy, as their production otherwise gives the strongest bounds. 
Another possibility is that the superpartners are not so heavy, but are difficult to detect because of 
a compressed mass spectrum, leading to much smaller visible energy than in mSUGRA benchmark 
cases. For our purposes here, compressed SUSY refers to the situation in which the mass ratio 
between the lightest super symmetric particle (LSP) and the gluino is significantly smaller than the 
prediction m^sp/rrig ~ 1/6 of mSUGRA. 

In a previous paper [2J, we investigated the reach of 2010 data from ATLAS, consisting of 35 
pb _1 of collisions at \fs = 7 TeV, for several classes of compressed SUSY models. In the present 
work, we will update this analysis to correspond to the 1.04 fb _1 data set analyzed by ATLAS in 
ref . [3j . Not only does this represent a huge increase in integrated luminosity, but also a revised set 
of signal regions compared to the 2010 data set analyses. The ATLAS analysis presents exclusion 
results for mSUGRA models, and for simplified models containing only squarks and a gluino but 
with the LSP mass held fixed at 0. In both cases the models tested are very far from the compressed 
case in which the mass difference between the gluino and the LSP is smaller. Our aims here are 
to see how the exclusions found for mSUGRA and simplified models translate into exclusions on 
compressed SUSY, for the various ATLAS signal regions, with particular attention to the exclusions 
that can be made in the most difficult case of very high compression. Other recent studies that 
include compressed SUSY and other non-mSUGRA searches at the LHC in a similar spirit can be 
found in 



The rest of this paper is organized as follows. Section [TT] describes the signal regions and our 
procedures. Section ILTT1 describes four classes of models, each of which depends on two parameters, 
the gluino mass M g and a compression parameter c, which are independently and continuously 
dialed to vary the overall superpartner mass scale and the ratio of gluino to LSP masses. Section 
IIVI gives results for the acceptances for these models with the various signal regions, and estimated 
exclusions based on the 1.04 fb _1 data set. Section |V] contains some concluding remarks. 



II. PROCEDURES AND SIGNAL REQUIREMENTS 



For this paper, we used the same tools as in our earlier work MadGraph/MadEvent 4.4.62 
23j was used to generate hard scattering events using CTEQ6L1 24| parton distribution func- 



tions, Pythia 6.422 25] for decays and showering and hadronization, and PGS 4 [26] for detector 
simulation. In SUSY models with compressed mass spectra, it is important to correctly gener- 
ate jets beyond the hard scattering event, by matching correctly (without overcounting) between 
matrix-element and showering/hadronization software generation of additional jets. We did this 
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> 130 
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Number of jets n 


> 2 


> 3 


> 4 


> 4 


> 4 


PT (jn) [GeV] 


> 40 


> 40 


> 40 


> 40 


> 80 


m cS [GeV] 


> 1000 


> 1000 


> 500 


> 1000 


> 1100 


Ef ss /m efi 


> 0.3 


> 0.25 


> 0.25 


> 0.25 


> 0.2 


ATLAS a x Acc [fb] 


< 22 


< 25 


< 429 


< 27 


< 17 



TABLE I: Summary of cuts for the signals A, B, C, D, E simulated here, following the ATLAS 2011 data 
analyses for 1.04 fb -1 Also shown on the last line are the ATLAS 95% CL bounds on the non-Standard 
Model contribution to the cross-section times acceptance in the five signal regions. (In the case of signal 
region E only, the m c ff > 1100 GeV requirement involves a sum over all jets with > 40 GeV, but the 
m e ff used in the E™ lss /m e g cut is a sum over only the leading 4 jets.) 



by generating each lowest-order process together with the same process with one additional jet at 
the matrix-element level, followed by MLM matching with P^-ordered showers with the shower- Kj< 
scheme with Q cu t = 100 GeV, as described in [27J and implemented in the MadGraph/MadEvent 
package. (It is much more time-consuming to include up to two extra jets at the matrix-element 
level. We found with some sample testing that it did not make a significant difference with our 
setup even for very compressed superpartner mass spectra.) For the detector simulation, we used 
the default ATLAS-like parameter card file provided with the PGS distribution, but with a jet cone 
size of AR = 0.4. Cross-sections were normalized to the next-to-leading order output of Prospino 

2.1 H- 

To define signals, we follow (a slightly simplified version of) the ATLAS cuts for multijets+i?™ lss 
from ref. [3j] • The signal requirements are summarized in Table HI Events are required to have at 
least one jet with pr > 130 GeV. The signal regions A, B, C, and D also require at least 2, 3, 
4, and 4 jets with pt > 40 GeV, respectively, while signal region E requires at least 4 jets with 
Pt > 80 GeV. These jets must have |r/| < 2.5. The leading three jets are required to be isolated 
from the missing transverse momentum according to A(fi(p™ lss , j) > 0.4. The effective mass m e g 
is defined as the scalar sum of the E™ 1SS and the pt's of: the leading 2 jets for A; the leading 3 
jets for B, the leading 4 jets for C, D; and all jets with pj- > 40 GeV for E. Then m e g is required 
to exceed 1000, 1000, 500, 1000, and 1100 GeV respectively for signal regions A, B, C, D, E. In 
addition, a cut is imposed on the ratio E^ ms /m e s of 0.3, 0.25, 0.25, 0.25, and 0.2 for A, B, C, D, 
E respectively. (For signal region E, only the 4 leading jet pt's are summed over in the m e s used 
in the E™ 1SS /m e ^ cut, even though the m e g cut uses an inclusive sum over jets.) Note that these 
cuts automatically imply a lower limit on E™ ss of 300, 250, 125, 250 GeV for signals A, B, C, D, 
respectively. For signal region E, a cut E™ 1SS > 130 GeV is imposed, although on an event-by- 
event basis this is usually superseded by the E^ lss /m c g cut. There is a veto of events with leptons 
I = (e,fi) with \rj\ < 2.4 and (2.47) for muons (electrons), and p^ T > 20 GeV that are farther than 
AR = a/ (Arj) 2 + (A</>) 2 > 0.4 from the nearest jet. Also shown on the last line of Table [J are the 
ATLAS 95% CL limits on non-Standard-Model contributions to the signal regions after acceptance 



and efficiency, as reported in ref. [3|. These will be used below to estimate the reach for compressed 
SUSY models. ATLAS also has searches requiring higher jet multiplicities 29J, leptons [30], and 
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b tagging [3l( , but these searches give significantly less reach for the compressed SUSY models we 
consider. Comparable searches by CMS have been reported in [32], [33j, and [34]]; we choose to use 
the ATLAS results only for reasons of convenience and familiarity. 

Because our tools for generating SUSY signal events and simulating detector response are not 
the same as those used by ATLAS, the cross-section and acceptance results found below clearly 
cannot be interpreted in exact correspondence to the ATLAS ones. However, we have checked 
that the results of our analysis methods correlate well to those in ref. [3j for a sample of mSUGRA 
models used there. For mSUGRA models with tan j3 = 10, Aq = 0, mo = 100 GeV, and m\ii = 180, 
210, 240, 270, 300, 330, 360, 390, 420, 450, and 480 GeV, we found agreement with the ATLAS 
acceptance*efficiency to be typically better than 15%, while for the same parameters but tjiq = 660 
GeV, the agreement was usually at the 30% level or better. Keeping these inevitable differences 
in mind, at least an approximate estimate of the true detector response may still be gleaned from 
the results below, and the general trends should be robust. 



III. COMPRESSED SUSY MODELS 

In this section, we define the compressed SUSY models for our study. Following our earlier 
work ref. Q] , let us parameterize the electroweak gaugino masses at the TeV scale in terms of the 
gluino physical mass as 



l + 5c\ „, A + 2c 
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Mi = ( — ^— ) M~ g , M 2 = — — M- g . (3.1) 



Here c parameterizes the degree of compression. The value c = gives an mSUGRA-like mass 
spectrum with gaugino masses approximately equal at Mqut = 2.5 x 10 16 GeV, and c = 1 gives a 
completely compressed spectrum in which the gluino, wino, and bino masses are equal at the TeV 
scale. The gluino mass M g is treated as the independent variable input parameter that sets the 
superpartner mass scale. We also select tan /3 = 10 and positive \i = M g + 200 GeV to compute 
the physical masses of charginos Cj and neutralinos Aj. For our first class of models, we take the 
first- and second-family squark masses to be: 

niu R = mg R = niu L = 0.96Mp, m| = m\ L - cos (2/3) m^, (3.2) 

and sleptons are taken degenerate with the squarks (so too heavy to appear in chargino and 
neutralino two-body decays). The top squark masses are taken to be = M g — 160 + c(180 — 
0.09M g ) and m ( - 2 = M g + 25, in GeV. The lightest Higgs mass is fixed at m h o = 115 GeV, and 
the heavier Higgs masses with m^o = 0.96M^. These choices provide relatively smoothly varying 
branching ratios as the compression parameter c is varied, although transitions of N 2 and C\ decays 
from on-shell to off-shell weak bosons are inevitable as the compression increases. In particular, 
the reason for the choice for the parameterization of the stop masses is to avoid suddenly turning 
on or off any 2-body decay modes as the parameter c is varied within each model line, by making 
sure that the the gluinos cannot decay to stops by kinematics for any of these models. The choices 
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FIG. 1: The masses of the 
most relevant superpartners 
for the models of class I de- 
fined in section ITTll as a func- 
tion of the compression pa- 
rameter c, for fixed M g = 
700 GeV. The case c = 
corresponds to an mSUGRA- 
like model. 
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c (compression parameter) 



for tan f3 and \x are arbitrary, and not very much would change if they were modified (within some 
reasonable range). We refer to the two-parameter class of models spanned by varying (M g ,c) as 
models I. 

The masses of the most relevant superpartners are shown in Figure Q] for the case M g = 700 
GeV, illustrating the effect of the compression parameter c on the spectrum. In this class of models, 
gluino and squark production dominate at the LHC. The gluino decays mostly by the two-body 
mode g — > qq or qq, and right-handed squarks decay mostly directly to the LSP, q~R — > qN\, while 
left-handed squarks decay mostly to wino-like charginos and neutralinos, qi — > q'Ci and qN2- 
The latter decay through on-shell or off-shell weak bosons: C\ — >■ W^N\ and iV 2 — > Z^N±, 
or N2 — > hNi when it is kinematically allowed. The visible energy in each event from these 
decays clearly decreases as the compression factor c increases, because of the reduction in available 
kinematic phase space. 

We define a second class of "heavy squark" models, II, which are the same as above but with 
all squarks taken very heavy, Mq = Mg + 1000 GeV. Thus when c = 0, the model classes I and 
II correspond approximately to mSUGRA with small and large mo, respectively. In these heavy 
squark models, the most important production cross-section is from gluino pair production, with 
subsequent gluino decays g — > C\qq' and N^qq and N±qq, with the first two typically dominating. 
The wino-like states then decay through on-shell or off-shell weak bosons, depending on the mass 
difference from the compression: Ci -> WWjVi and N 2 -» Z^N X or hNi, with the last dominating 
if kinematically allowed. 

One motivation for compressed SUSY is that taking M3/M2 much less than the mSUGRA 
prediction can significantly ameliorate 0, 35, 36] the supersymmetric little hierarchy problem. 
This provides motivation to consider compressed SUSY models in which winos are heavier than the 
gluino. Therefore, we define a third and fourth class of models, III ( "heavy winos" ) and IV ( "heavy 
winos and squarks"), to be the same as models I and II respectively, but with M2 = M g + 100 GeV 
in each case. In models III, all first- and second-family squarks decay directly to the LSP: q — > qN\, 
while in models IV the squarks decouple from the discovery or best exclusion limit processes. The 
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gluino has direct two-body decays to quarks and squarks as before. 

The model classes I, II, and III were exactly those we used in ref. |2| in the context of limits 
obtainable with 35 pb _1 at LHC, while the model class IV corresponds approximately to the heavy 
squark limit of the simplified gluino/squark models in Q], but with neutralino LSP masses that 
are here non-zero and vary continuously with c. Thus the models that we discuss here provide a 
quite different slicing through the MSSM parameter space than those found in the experimental 
collaboration papers. We now proceed to use them to examine how the ATLAS exclusions on 
non-Standard Model cross-section times acceptance times efficiency impact the parameter space 
for SUSY with compressed mass spectra. 



IV. RESULTS OF SIMULATIONS 

In compressed SUSY models, the visible energy in jets is reduced compared to mSUGRA mod- 
els, leading to a significant reduction in acceptance for signal events. In Figure [2j we show the 
fractional acceptances after all cuts for the models in class I with c = — 0.1, 0, 0.1, . . . , 0.9 and 
Mg = 300, 400, . . . , 1300 GeV. Results are shown for each of the signal regions A, B, D, and E. (We 
find that signal region C is not competitive for setting limits in any of the models we consider, see 
Figure H] below, so it is not included in Figure [2j) 

For fixed values of the mass difference M g — Mlsp, the acceptance tends to approach a limit 
for sufficiently large Mg. Conversely, for sufficiently large M g — Mlsp, the acceptance tends to be 
relatively flat, but falls dramatically for M g — Mlsp ^ 450 GeV for signal regions A, B, and D, 
and for an even larger range of the mass difference for the high-mass signal region E. For severe 
compression M g — Mlsp < 150 GeV, the signal regions A and B can be seen to retain acceptance 
more than the signal regions D and E do, although in each case the acceptance declines to well 
below 1% for the most extreme compression, even when the gluino mass is very large. 

The acceptances for the heavy squark models of class II are similarly shown in Figure [3] for 
gluino masses from 300 up to 1000 GeV. Here the acceptance tends to increase more steadily with 
increasing Mg. There is sometimes a notable decrease in acceptance for signal regions A and B 
as the mass difference M g — Mlsp increases, so that the largest acceptances are achieved with 
non-zero compression, that is, when M g — Mlsp is n °t maximum. This perhaps surprising effect, 
studied in [jj, is due to the fact that as the compression increases, the m e g distribution becomes 
soft faster than the i?™ ss distribution does, so that more events pass the E^ lss /m e s cuts. 

The acceptances for models in classes III (heavy winos) and IV (heavy squarks and winos) are 
qualitatively similar, and so are not shown. 

In Figure HJ we show contours of cross-section times acceptance for all five signal regions, for 
each of the model classes I, II, III, IV in separate panels. The contours for each signal region are 
for the corresponding ATLAS limits on non-Standard Model processes listed in Table U (taken from 
0]), so that the regions to the left of each contour may be regarded as the approximate exclusion 
regions for that signal definition. In the panels for model classes I and II, the (orange) dotted line 
indicates the case c = in which the ratio of gaugino masses at the TeV scale is approximately the 
same as mSUGRA. For the light squark class of models I, the best signal regions for exclusion are 
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FIG. 2: Acceptances for the four signal regions A, B, D, E defined in Table HI for the class of models I 
defined in the text. The lines on each graph correspond to different values of the gluino mass. The dots on 
each line correspond to compression factors c = — 0.1, 0,0. 1, ...0.9 from right to left. 



A (when M g — Mlsp > 400 GeV) and B (for most of the range of smaller mass differences, except 
for the most extremely compressed case). 

Signal regions A and B likewise give the best exclusions for models in class III (heavy winos but 
light squarks). Because the gluino and squark decays in this class of models do not pass through 
the intermediate cascade step of winos, the visible energy per jet tends to be larger, leading to 
stronger exclusions as shown. In both of the model classes I and III with squarks slightly lighter 
than the gluino, we find that even in the case of extreme compression one can still set a limit of 
better than M g > 600 GeV. 

In models with heavy squarks, the limits are much worse, since the main SUSY production is 
only gluino pairs. For models of class II with heavy squarks and light winos, the best limit is set 
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FIG. 3: As in Figure [2j but for the models II (heavy squarks) defined in the text. 



using the high-mass signal region E when M g — Mlsp > 200 GeV. For smaller mass differences, the 
3-jet signal region B sets better limits, because the m e s distribution for the signal events becomes 
too soft. Still, it should be possible to set a limit of about M g > 450 GeV using signal regions 
A and B, even in the case of extreme compression with c = 0.9. Qualitatively similar statements 
apply to models in class IV with both squarks and winos decoupled. Note than in all cases, signal 
region C is comparatively ineffective in setting limits, because the backgrounds are too large. 
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FIG. 4: Contours of constant cross-section times acceptance for the five signal regions defined in Table HI 
in the Mg — Mfj vs. Mg plane obtained by varying the gaugino mass compression parameter c between 
—0.1 and 0.9. The four panels correspond to the models I (light squarks), II (heavy squarks), III (heavy 
winos), and IV (heavy winos and heavy squarks). The dotted lines in the first two cases corresponds to the 
mSUGRA-likc case c = 0. Increased compression is lower in each plane. 
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FIG. 5: The impact of systematic uncertainties in the signal rates on contours of constant cross-section 
times acceptance for selected signal regions. The middle line in each case is the same result as in Figure 
SJ The corresponding left and right contours show the impact of decreasing and increasing the total signal 
production rate by 25%. The two panels correspond to the models I (light squarks, left) and II (heavy 
squarks, right). 



The preceding results were all obtained using the Prospino NLO default renormalization and 
factorization scale choices and without taking into account the possible systematic uncertainties in 
the signal production cross sections. In general, uncertainties in QCD production cross-sections are 
notoriously difficult to estimate. It is well-known that variation of renormalization and factorization 
scales and PDFs do not give reliable estimates of the production cross-section uncertainties. To 
illustrate the potential impacts of these uncertainties, we show in Figure [5] how the results vary 
when changing the assumed total signal cross-section by ±25% for the models of class I and II. Only 
the two signal regions that give the strongest limits over the most significant ranges of M g — Mlsp 
are shown in each case. This variation results in a change in the gluino mass limit in these models 
that can exceed ±50 GeV, depending on the model parameters. 



V. OUTLOOK 



In this paper, we have studied the reach of 1 fb _1 of LHC data at yfs = 7 TeV for compressed 
SUSY models, extending our earlier results for 35 pb _1 in [2]. We found that even in the most 
compressed case studied, in which the gluino is only about 9% heavier than the LSP, the limit on the 
gluino mass should be about M g > 600 GeV for squarks that are slightly lighter than the gluino, and 
about Mg > 450 GeV when squarks are very heavy. The best limits (and discovery potential) come 
from signal regions which require 2 or 3 jets. In designing future searches for compressed SUSY, 
it is probable that the best reach will be obtained by increasing the cut on E^p lss as necessary 
to reduce the backgrounds, rather than by very hard cuts on m e g- (or Ht)- This is because as 
the compression increases, both the i?™ lss and m e s distributions get softer, but the latter more 
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drastically. (A more precise quantitative statement about this is beyond the scope of this paper, 
since it would require detailed background estimates including crucially detector response-specific 
information.) Future searches should take into account that signal regions optimized for mSUGRA 
and for simplified models with massless or light LSPs will therefore not do very well for compressed 
SUSY models, and this effect will become more significant as higher mass scales are probed. 
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